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Abstract — In cognitive radio systems, cyclostationary feature 
detection plays an important role in spectrum sensing, especially 
■ in low SNR cases. To configure the detection threshold under a 
I certain noise level and a pre-set miss detection probability Pf, it's 
, important to derive the theoretical distribution of the observation 
variable. In this paper, noise distribution in cyclic frequency 
domain has been studied and Generalized Extreme Value (GEV) 
I distribution is found to be a precise match. Maximum likelihood 
estimation is applied to estimate the parameters of GEV. Monte 
Carlo simulation has been carried out to show that the simulated 
ROC curve is coincided with the theoretical ROC curve, which 
proves the efficiency of the theoretical distribution model. 

Index Terms — Cognitive radio, spectrum sensing, cyclic feature 
detection, noise distribution on cyclic frequency domain 

! L Introduction 

Spectrum sensing plays an important role in cognitive radio 
(CR) systems, and cyclostationary feature detection is one 
, of the main technologies for spectrum sensing in low SNR 
' cases [1-2]. The idea of cyclic detection is that one CR node 
samples the RF signals, transform the time domain signals 
into cyclic frequency domain, then decide primary user's 
', occupancy of a target band regarding whether the cyclic 
' spectrum on a significant cyclic frequency is above a certain 
I threshold [2]. 

Although the simulation results of cyclostationary feature 
detection has been studied [3-5], less attention has been 
paid to analyze the noise distribution on cyclic frequency 
domain. Consequently, theoretical function between detection 
threshold \th and miss detection probability Pf is not 
available, which makes it difficult for practical system design. 
For comparison, the theoretical function between detection 
threshold Xth and miss detection probability Pf of energy 
detection is available by using central and non-central chi- 
square distribution to model the distribution of the observation 
variable affected by time domain Gaussian noise [6], which 
makes energy detection to be a practical method for spectrum 
sensing. However, its performance in low SNR cases is much 
poorer than cyclostationary feature detection [7]. In this 
paper, noise distribution on cyclic frequency domain has been 
analyzed and generalized extreme value distribution is found 
to be a precise match of the observation variable affected 



by time domain Gaussian noise. A fast cyclic frequency 
domain feature detection algorithm [7] has been introduced 
to evaluate the coincidence between theoretical ROC curve 
and the simulated ROC curve, which proves the reliability 
of theoretical distribution model and feasibility of practical 
system design. 

The rest part of the paper is organized as follows: Section 
II describes the system model of cyclostationary feature 
detection. Noise distribution on cyclic frequency domain 
is analyzed in section III. A fast cyclic frequency domain 
feature detection algorithm has been introduced in section IV. 
Simulation results are given in Section V. Finally, conclusions 
are drawn in Section VI. 



II. SYSTEM MODEL 

The spectrum sensing problem can be modeled as hypoth- 
esis testing. It is equivalent to distinguishing between the 
following two hypotheses: 

r Ho : yit) = n{t) 

\ Hi : y{t) = x{t) + n{t) ^ ' 

y{t), x{t) and n{t) denote the received signal, the primary 
user's transmit signal and the Gaussian noise, respectively. Hi 
and Hq represent the hypothesis that the primary user is active 
or inactive. Due to the existence of noise, a certain threshold 
Xth should be set to decide whether a primary user is active 
or not. Probability of detection (Pd) and false alarm (Pf) are 
defined to evaluate the detection performance: 

/ Pd=p{y, > Xth\Hi) 

\ Pf = piy, > \th\Ho) 

The goal of detection is to maximize Pd while maintain a 
given Pf. 

When feature detection is applied, the detection model (1) 
changes into: 

J Ho ■■ S^{f) = S'"(/) 

\ Hi:S^if) = S^if) + S^,if) ^ ' 



Sy{f) is the spectrum correlation density (SCD) of the re- when k — 
ceived signal y{t), S'^{f) and S^^{f) is the SCD of x{t) and 
n{t), respectively [8]. 

Theoretically, Gaussian noise n{t) is not a cyclostationary 
statistic process, then S"{f) = when a 7^ [8]. As for 
cyclostationary signal x{t), there is a significant frequency 
set {ao}, on which S""{f) ^ 0. Due to the ideal non- 
cyclostationary characteristic of Gaussian noise, any pre-set 
threshold \th on a significant frequency will lead to = 1 



and Pj = 0. 



III. Noise distribution on cyclic frequency 

DOMAIN 

In practice, SCD is calculated for limited length signals, 
therefore, S"{f) ^ when a 7^ [9]. As shown in Fig.l, 
the background noise is obvious on / ^ a square when 
calculating SCD of a noise interfered AM modulated signal. 




Cyclic Fraquencrt^yHz 



Fig. 1. Cyclic Spectmrn Density of AM at SNR=-10db 

In order to analyze the background noise distribution, a 
limited-length Gaussian noise sequence n*^(j) is considered: 



n{iK + j),0 < j < K -1 
0, else 



(4) 



where K is the length of the analysis window, z = 0, 1, 00 
is the index of analysis window. Noise data in each analysis 
window are transformed into cyclic frequency domain, and 
then mapped from f ^ a square to a axis through the 
following expression: 



7V*^ (a) 



max|5,'^(/)| 



(5) 



For each cycle frequency ao, the cyclic spectrum value is 
aligned to a set {Ni = N"^^ (ao)}, z = 0, 1, L. According 
to the extreme definition of N'^^ (ao) in (5), Generalized 
Extreme Value (GEV) distribution is adopted to model the 
cyclic frequency domain noise [10]. The density function is: 
when K 7^ 



/(a;|K,^,cr) 



4)= 



'(1 



(x - n) 



f{x\0, ^, a) = {-)exp{--exp{-- — -) - — -)) (7) 
a a a 

where k is the shape parameter, /i is the position parameter, 

cr > is the scale parameter. And the parameters k , /i , cr 

can be estimated by maximum likelihood estimation based on 

the noise sequence {Ni}. For most cases, k « 0, then the 

likelihood function is defined as follows: 

/(^,a) = -LZna-V(^^)-Vexp{-(^^)} (8) 
Let # = 0,# = then: 

Ofl ' OCT 



Eti(^.-M)(l-e-(^-«/*) = i<7 



(9) 



where fi and a are the estimated value of /i and a. By solving 
(9), jl and a are obtained. 

After that, the likelihood function for k is defined as: 



1 ^ 

1{k, fi, a) = -Nina - (1 + -) X! ^""[^ + 



N 



n.; — U^^ 1 



)] 



(10) 



Let ^ = 0, the estimation of k can be obtained by solving 
(10). " 

IV. Fast Cyclic Frequency Domain Feature 
Detection Algorithm 

Fig. 2 shows the block diagram of a cyclic frequency domain 
feature detector [7]: 



SCD Calculation 



Cyclic Frequency 
Component 



Fig. 2. Block diagram of a feature detector 

In this detector, time domain signals are transformed to 
cyclic frequency domain, and then mapped to a axis and 
get extreme value for each a, finally, compared with a pre- 
set threshold \th to determine the occupancy of primary 
user. To decrease the computational complexity of feature 
detection, only one cycle frequency ao for a modulated signal 
is calculated for spectrum sensing [7]. The probability of 
detection iPd) and false alarm (P/) is defined by: 

Pd = PriT > \th\Hx) 
Pf = Pr{T > Xth\Ho) 

where Pr( ) is CDF of generalized extreme value distribution. 
Substitution of (6) and (7) into (11) yields: 



(11) 



(6) 



Pf ^ PriT > Xth\Ho) 



f{x\k,jl,a)dx (12) 



For a pre-set Pf, the threshold can be estimated as: 

A.. = ('?-^,^^-^?-^)'''^' (13) 
y^jL- crlogyp,K = 

where yp = — /og(l — Pf). 

V. Simulation Results 

In this section, Monte-Carlo simulation results are presented 
to prove the reliability of the upper analytical results between 
Pf and the threshold Xth- 

Frequency smoothing method in [11-12] is applied to 
estimate SCD of a time domain noise signal. Simulation 
parameters are listed in TABLE 1. Length of the analysis 
window is set to be 4096 and totally 10000 cyclic spectrum 
values on cyclic frequency ao = 2fc, aligned by analysis 
window index, are considered to evaluate the theoretical curve. 

The histogram of the aligned cyclic spectrum values is 

TABLE I 
Simulation parameters list 



Parameter 


Value 


Modulation type 


AM 


Carrier frequency 


1 MHz 


Bandwidth 


10 KHz 


Sampling frequency 


3 MHz 


Sampling time 


1.365ms 


Channel 


AWGN 


Window type 


hamming 


Frequency smoothing length 


1300 


Sampled data length(K) 


4096 



shown in Fig 3, with compared to GEV distribution. It is 
proved that the GEV distribution precisely match the cyclic 
frequency domain noise data. 



receiver operating characteristics (ROC) , which are theoret- 
ically derived from GEV distribution, are plotted to compare 
with those derived from Monte-Carlo simulation. As to the 
theoretical curve, Pf is pre-set according to system require- 
ment. By using (13), a theoretical threshold Af^ is obtained. 
Finally, received signals under hypothesis Hi are compared 
with the threshold to obtain the statistics result of Pd- As to 
the simulated curve, threshold are chosen to be the same as the 
theoretical threshold {At/i}, after that, received signals under 
hypothesis Hi and iJo are compared with each Af^ to obtain 
the statistics results of Pd and Pf. Finally, plot these {Pd,Pf) 
points to form a continuous curve. 

Experiments results are shown in Fig.4, we can see that the 
theoretical ROC curve (red highlighted) precisely match the 
simulated ROC curve (green highlighted) for different received 
signal power levels. It is proved that the generalized extreme 
value distribution is efficient to model the noise distribution 
on cyclic frequency domain. 
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Fig. 4. ROC curves under different signal power levels 
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Fig. 3. GEV curve for the noise samples 

For further proof of the proposed model, the curves of 



VI. Conclusion 

In this paper, noise distribution on cyclic frequency domain 
is studied and generalized extreme value (GEV) distribution 
is found to be an efficient method to model the cyclic 
frequency domain noise. Maximum likelihood estimation is 
applied to estimate the parameters of GEV. Sensing threshold 
is consequently derived from system requirements (a pre-set 
Pf) and theoretical CDF of GEV distribution. Monte Carlo 
simulation has been carried out to prove that the simulated 
ROC curve is precisely coincided with the theoretical ROC 
curve. 
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